Glass Formation in a Periodic Long-Range Josephson Array 
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We present an analytic study of a dynamical instability in a periodic long-range Josephson array 
frustrated by a weak transverse field. This glass transition is characterized by a diverging relaxation 
time and a jump in the Edwards- Anderson order parameter; it is not accompanied by a coinciding 
static transition. 
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Glass formation in the absence of intrinsic disorder is a 
long-standing problem. Although vitrification is ubiqui- 
tous, a minimalist microscopic model of this phenomenon 
remains the subject of active discussion. jjj Because glass 
formation is a dynamical transition that is not necessarily 
accompanied by a static one, it lies outside the framework 
of the Landau theory. Furthermore this glass transition 
leads to a low-temperature phase with broken ergodicity 
without the selection of a unique state. A successful phc- 
nomenological theory should describe the "partitioning" 
of phase space below the transition temperature (Tq) into 
an exponential number of metastable states; specifically 
it should explain how the system becomes "stuck" at Tq 
in one of these states that is separated from the others 
by barriers that scale with the system size. 

Unfortunately a basic theory of glass formation has not 
yet been found. Several microscopic non-random models 
have been proposed; most were studied via a mapping to 
disordered systems . Recently possible glassiness in the 
absence of disorder has been studied in a periodic long- 
range Josephson array using a direct, analytic approach 
||; furthermore this system may be realized experimen- 
tally. An analysis of its static behavior indicates a first- 
order transition into a low-temperature phase character- 
ized by an extensive number of states separated by infi- 
nite barriers. In this Letter we continue the study of this 
system and show that it displays a true dynamical insta- 
bility that precedes the static transition, as expected in 
a glass. 

The proposed array is a stack of two mutually perpen- 
dicular sets of N parallel wires with Josephson junctions 
at each node (Figure 1) that is placed in an external 
tranverse field. The classical thermodynamic variables 
of this system are the superconducting phases associated 
with each wire. Here we shall assume that the Josephson 
couplings are sufficiently small so that the induced fields 
are negligible in comparison with H. We can therefore 
describe the array by the Hamiltonian 
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with Jot. = 



-exp(2majk/N) and 1 < (j, k) < N 



where j(k) is the index of the horizontal (vertical) wires; 
s m = e z ^ m where the (j) m are the superconducting phases 
of the wires. Here we have introduced the flux per 
unit strip, a — N HI 2 / '<po, where I is the inter-node spac- 
ing and (f>o is the flux quantum; the normalization has 
been chosen so that Tq does not scale with N. 




where J mn is the coupling matrix 



Fig 1. The phase diagram of the array (inset) where 
Tq indicates the temperature associated with the dynam- 
ical instability discussed in this Letter, T$ is the specu- 
lated equilibrium transition temperature and T m is the 
"superheating" temperature where the low-temperature 
metastable states cease to exist. 

Because every horizontal (vertical) wire is linked to 
every vertical (horizontal) wire, the number of nearest 
neighbors (z) in this model is N; we can therefore study 
it with a mean-field approach. For -k <C ot < 1 the 
number of low-temperature metastable solutions is ex- 
tensive p| . Furthermore this degeneracy develops simul- 
taneously with the instability of the paramagnetic phase; 
at this temperature interactions do not favor a particular 
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state and, because z ~ N, the barriers separating these 
low-temperature metastable solutions are effectively infi- 
nite. The static response displays no soft mode instabil- 
ity but indicates a jump in the Edwards- Anderson order 
parameter in the vicinity of Tq = -^S= . 

Before presenting our quantitative treatment of the 
dynamical behavior of this array, we discuss the quali- 
tative picture of the glass transition that emerges from 
our results (cf. Fig. f). As T approaches T+, where 
T+ ~ To, there appear a number of metastable states in 
addition to the paramagnetic free-energy minimum; most 
likely they are energetically unfavorable and thus do not 
"trap" the system upon cooling from high temperatures. 
As T -> T+, the paramagnetic minimum is "subdivided" 
into an extensive number of degenerate metastable states 
separated by effectively infinite barriers, and the system 
is dynamically localized into one of them. Qualitatively, 
in the interval T rn > T > Tq there appear many local 
minima in the vicinity of the paramagnetic state sepa- 
rated by finite barriers; these barriers increase contin- 
uously and become infinite at T = Tq. Each of these 
minima is characterized by a finite "site magnetization" 
mi = (si/T where"site" refers to a wire. When T > Tq 
thermal fluctuations average over many states so that 
(rrii) = 0. At T = Tq the system is localized in one 
metastable state and there is an associated jump in the 
Edwards- Anderson order parameter, (q — jj Si( TO i) 2 )- 
The low-temperature phase is characterized by a finite q 
and by the presence of a memory, lirat'_>oo A(t, f') 7^ 
where A(t, t 1 ) is the anomalous response. We expect that 
at T = Tq, the metastable states are degenerate and thus 
there can be no thermodynamic selection. However at 
lower temperatures interactions will probably break this 
degeneracy and select a subset of this manifold; then we 
expect an (t — ► 00) equilibrium first-order transition (Ts) 
which should be accompanied by a jump in the local mag- 
netization. In order to observe this transition at Ts the 
array must be equilibrated on a time-scale (tw) longer 
than that (t^) necessary to overcome the barriers sepa- 
rating its metastable states; tA scales exponentially with 
the number of wires in the array. Thus the equilibrium 
transition at Tg is observable only if tw — ► 00 before the 
thermodynamic limit (N — > 00) is taken; in the opposite 
order of limits only the dynamical transition occurs. 

We now begin a more quantitative analysis of the dy- 
namic instability in this periodic array. Because our fo- 
cus is on the long-time behavior of this system, we ex- 
pect the details of the single-spin dynamics to be irrel- 
evant; we therefore choose to study the simplest form, 
namely that of soft spins with Langevin relaxational dy- 
namics. More specifically, we introduce a "potential" , 
V(Si) = Vq(\S\ 2 — l) 2 , at each wire which constrains the 
magnitude of each spin, \Si\ ~ 1, and assume the equa- 
tions of motion 
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where tj, is a microscopic time-scale. The dynamics (|3|) 
reproduces the dynamics of the overdamped Josephson 
junctions with individual resistance R if 77, = w^bT 



and Vq — > 00. In order to average the solution of (|^) over 
the thermal noise, £, we use a generating functional. For 
example, the average supercurrent in the array is given 
by 



(Iij) = J IijexpA[S,S]T>SVS. 



(•5) 



Here the current 1^ = (§§) ImS* JijSj and the action is 
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where we have not included the terms that arise from the 
Jacobian since they do not affect the long-time response. 

We perform our calculations by resumming the terms 
in the -tt expansion of (§) which are leading order in -k; 
this is a dynamical analogue of the high-temperature se- 
ries expansion previously used to study the static behav- 
ior of this array. || The crucial ingredients of this tech- 
nique are the response (G mn (t,t') = (s m (t)s n (t'))) and 
the correlation (C mn (t,t') = (sm(t)s^(t')}) functions. 
For T > Tq, these depend only on the time-differences 
and thus can be related by the Fluctuation-Dissipation 
Theorem 



G <J -(t-t') = - 
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The leading diagrams (in -k) for Gij(t — t') are displayed 
in Fig. 2a. The presence of the "constraining poten- 
tial" V in the action (^) results in finite higher-order ir- 
reducible single-site spin correlations, which play the role 
of interaction vertices in this diagrammatic technique. 
However the corrections to the response function shown 
in Fig. 2b are small in in comparison with those in 
Fig. 2a. 

We emphasize that, as in the static case, the single- 
site response is renormalized; here we consider the local 
Green's function (G(t — t')) that is irreducible with re- 
spect to the Jij lines. Possible self-energy corrections to 
G(t — t') are shown in Fig. 2c and will be discussed be- 
low. Summing the geometric series shown in Fig. 2a, we 
obtain 
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for the response function connecting wires of the same 
type (horizontal/vertical). The matrix (J'J)ij depends 
only on the "distance" i — j and acquires a simple form 
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in Fourier space (jt J) p — (Jg /a)8(air — \p\); in this rep- 
resentation the Green function becomes 



gup) = - 



\p\) , 0(|p|-mr) 
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the high-temperature single-site nonlinear susceptibility 
X 3 = -^r. The set of equations, (§,§), © and @, 
are sufficient to determine the response and the correla- 
tion functions of the array. 



The static limit (u — 0) of TG" 1 coincides with the lo- 
cator, A(T), discussed previously P); in the absence of 
Onsager feedback terms Gq~ 1 = 1. Therefore we see in (jj^) 
that there would be a static instability at Gq 1 = G c = 
For 6 = (T - T ) /T > y/a all feedback effects arc 

negligible. Here the time-dependence of Gj 1 is set by a 
microscopic time-scale Tb and G^ (a>) = f3A(T) — iwrfe; 
inserting this G^ 1 (tj) into (||) we see that the long-time 
behavior of G w (p) is dominated by the first term which 
results in a long relaxation time r = (2/a)rb where 

j2 

a = (3 (A — » 28. In this regime the single-site 
response function is 



G(t) 



-t/r 
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t > T b . 



(10) 



1 so that it ap- 
■ and instead 



At lower temperatures, < y^a, the feedback ef- 
fects become important; they modify Gq 
proachcs G c only asymptotically at T — 
a first-order transition occurs ||. The retardation of 
the Onsager terms is also crucial and significantly af- 
fects the long-time behavior. Qualitatively the resulting 
dynamical instability, described below, is due to the time- 
dependence of the cavity field which itself is determined 
by single-site susceptibilities; the time-scale associated 
with the relaxation of G(t) increases continuously due to 
feedback through the Onsager terms. Formally the latter 
introduce an additional frequency-dependent part of the 
local response 
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as a self-energy E w such that Eo = since we have chosen 
our normalization so that Gq 1 = f3A(T). 

The on-site self-energy terms (see Fig. 2c) are the sim- 
plest for a <C 1, and thus we will consider this regime. 
We focus on the long-time response of this sytem which 
is dominated by the first term in (|9|) when Gq 1 « G c ; its 
weight is proportional to a (cf. ([l0|)). Thus in the limit 
of a <C 1, the slowly decaying parts of D(t) and G(t) 
scale with a and the dominant self-energy contribution 
contains the minimal number of these functions (X^ 3 -* in 
Fig. 2c) and is given by 



D 2 (t)G(t) (e 1 ' 



l)dt 



(12) 



where D and G are single-site correlation and response 
functions respectively. Q Here T is the four-spin vertex; 
we neglect its transient time-dependence and approxi- 
mate it by its static value T — — 1 which is determined by 
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Fig 2. Diagramatic expansion for the response function 
G(t); the dashed, solid and thick (with dot) lines are the 
coupling matrix J, the single site irreducible (G(t)) and 
the full response (correlation) functions respectively, a. 
Leading order in 1/N. b. Subleading order in 1/N. c. 
Leading terms in the expansion for G(t); the first dia- 
gram dominates at small a. 

Since we would like to detect a dynamical instability, 
we only consider the long-time behavior of the response 
function, i.e. 



G(t) = a 



We make the ansatz 



2(E W + iuJT b ) 
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(14) 



and show that it solves the system of equations 
(@'l I! H)' m tne P rocess we a l so obtain the time-scales 
t and t\ controlling the physical response. We have 
chosen this form of E w so that Eo = and £(t) ~ 
exp(—t/ri). Inserting this ansatz into (|l3|), we obtain 
response and correlation functions that are simple expo- 
nentials for long times; therefore the resulting self-energy 
(cf. (p^)) is also exponential and thus of the same form 
as the initial ansatz. Equating the parameters of (|l|) 
and (f[2|), we find that the condition for self-consistency 
at a — > a c is 



T l = 77 



l n 



3 a — a c 



r = 2a c n, 



(15) 



The resulting long-time part of the single-site response 
and the correlation functions are 



G(i) = 



2a _^ 2a _^ 

-e t r , D(t) = — e t r 



3ar R 
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(16) 



3 



where tr — 3t%, and we see from (15 T^) that the longest 
physical time-scale, tr, diverges continuously at G = 
-|(2a) 1 / 4 (cf. Fig. 1). At 9 = 6 G the long-time part of 
D(t) shown in Jl6| ) becomes constant q — (2a) 1 / 4 show- 
ing that the Edwards- Anderson order parameter jumps 
at the Sg- The resulting phase diagram is displayed in 
Figure 1; we note that T G = (1 + 6 G )T occurs at a lower 
temperature than T m , where the last low-temperature 
metastable states disappear M as discussed above. 

The response function G(t) is a susceptibility with re- 
spect to the field conjugate to S = exp i(f>, and thus can- 
not be measured directly. However, using G(t) found 
above, we can determine the ac-response to a time- 
varying physical magnetic field H (t) which is experimen- 
tally accessible. We focus on the total magnetic moment 
of the array generated by the Josephson currents: 



M = \ (ID f Y,( Sm ^ mnS ^ 



(17) 



where ^J mn = imnj mn if m and n are indices referring 
to horizontal and vertical wires respectively. We would 
like to determine the response in this magnetization to a 
time-varying field; we use the fact that M. = for static 
H to write 



dM(t) 
dH(f) 



l 2 ReTrJG(t,t')JD(t',t). (18) 



In order to evaluate (|18| ) we will need the response 
function connecting wires of different type (horizon- 
tal/vertical) 



G= J^r- 



1 



G" 2 -/3 2 JtJ 
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and of the same type (cf. (||)). We use the Fourier repre- 
sentation of jt J and that of G and D to determine the 
ac-response, , in terms of the single-site response 

functions ( L = Nl): 
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We can insert the response and correlation functions 
found above to determine the ac susceptibility %uj — dA '' 
which leads to 



2N 
~~9~ 



a uj + 2i/tr 
where tr is the longest time-scale of the response: 
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(22) 
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From Eq. (|22|) we see that the divergent relaxation times 
are directly observable in the physical a.c. magnetic re- 
sponse of the array. The zero frequency limit of the a.c. 
susceptibility jumps to a finite value at T — Tq, indicat- 
ing the development of a finite superconducting stiffness 
at the transition. Therefore, the measurement of this a.c. 
response in a fabricated array serves as a probe of glass 
formation. 

In summary, we have presented a periodic model 
which displays a dynamical transition where the system 
"freezes" into one of an extensive number of metastable 
states. Due to their large degeneracy, these states will 
not be selected by a Boltzman weight under equilibrium 
conditions. This glass transition at Tq is characterized by 
diverging relaxation times and by an accompanying jump 
in the Edwards- Anderson order parameter; the phase di- 
agram of the array is displayed in Figure 1. It would 
be interesting to study the physical properties of this 
array below the glass transition temperature in its non- 
crgodic regime; in particular we expect "memory" effects 
in the form of an anomalous response function and "fin- 
gerprints" of the individual metastable states in its phys- 
ical behavior. Since any uncertainty in the position of the 
wires introduces randomness in this array, it also offers 
the opportunity to study the crossover between glasses 
with spontaneously-generated and quenched disorder. 
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